Final-State Interaction as the Origin of the Cronin Effect 
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Instead of adhering to the usual explanation of the Cronin effect in terms of the broadening of 
the parton transverse momentum in the initial state, we show that the enhancement of hadron 
production at moderate pr in d+Au collisions is due to the recombination of soft and shower 
partons in the final state. Such a mechanism can readily explain the decrease of the Cronin effect 
with increasing rapidity. Furthermore, the effect should be larger for protons than for pions. 



o 
o 



> 

O 
O 
m 
o 

o 



X 



PACS numbers: 25.75.-q,24.45.-z,24.85.+p 

The conventional explanation of the Cronin effect 0, 
i.e., the enhancement of hadron spectra at high pT in 
proton-nucleus (pA) collisions, is that it is due to multi- 
ple scattering of projectile partons by the target nucleus 
before the production of a minijet by a hard scattering 
Pi . All models on the effect are based on the traditional 
approach to hadron production at high px, which is to 
follow a hard-scattered parton by a fragmentation of that 
parton. Since in that paradigm there is nothing more one 
can do with the final state, all models focus on the initial 
state, and they differ only in the way the broadening of 
the intrinsic traverse momentum is implemented. In this 
paper we consider a drastically different approach to the 
problem. We do not assume initial-state broadening, and 
treat the hadronization process in the final state by re- 
combination It will be shown that the Cronin effect 
can be satisfactorily explained at all centralities. In so 
doing we eliminate the necessity of putting in by hand 
successive transverse kicks of appreciable magnitude in 
the initial state. 

The idea that final-state interaction may contribute to 
the Cronin effect is not new 0, . However, no theoret- 
ical model has ever been proposed to demonstrate that 
the idea can be translated into quantitative accounting 
of the effect. Here we work in the specific framework 
of the recombination model and make concrete predic- 
tions. Although the model has long been used to treat 
hadronization in the fragmentation region, a number of 
groups have recently found that recombination is more 
important than fragmentation at small and moderate px 
(0 < pt < 8 GeV/c) at midrapidity in heavy-ion colli- 
sions at ^/s = 200 GeV A central issue in the 
model has always been the determination of the appro- 
priate distributions of partons that recombine, and it has 
been the major problem to tackle for the application of 
the model at low pt in hadron-hadron |^ and proton- 
nucleus collisions. For high px we have very recently 
determined the distributions of shower partons created 
by a hard parton and then found that the recombi- 
nation of thermal and shower partons is crucial in repro- 
ducing the observed spectra at moderate px in Au+Au 
collisions It is the thermal-shower recombination 
that also holds the key to the Cronin effect in pA colli- 



sions. It is an effect due to the interaction between soft 
and hard partons during hadronization in the final state. 

Unlike the case of Au-|-Au collisions at RHIC where the 
hydrodynamical expansion of a hot, dense system that is 
created leads to a large body of thermal partons, one 
does not expect such a scenario in pA collisions. Never- 
theless, there are soft partons that take the place of the 
thermal partons. They can participate in the formation 
of hadrons at moderate px- The number of such soft 
partons decreases with increasing impact parameter and 
rapidity, so our explanation of the Cronin effect natu- 
rally leads to a reduction of the effect at lower centrality 
and hi ghe r rapidity, in accordance to the observation at 
RHIC 0. 

The inclusive distribution for the production of pions 
can be written in the recombination model, when mass 
effects are negligible, in the invariant form 
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where Fqq{pi,p2) is the joint distribution of a, q and q at 
Pi and p2 , and Rt^ (pi ,p2,p) is the recombination function 
for forming a pion at p: i?„(pi,p2,p) = {piP2/p)S{pi + 
P2 ~p)- Restricting p to the transverse plane, the dis- 
tribution dNT^/d'^pdy\y=o, averaged over all (/>, with px 
denoted by p, is |l3 | 
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dpiFqq{pi,p- pi). 
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This equation is applicable to any of the pp, pA and 
AB collision types; only Fqq depends on the colliding 
hadron/nuclei. In general, Fqq has four contributing com- 
ponents represented schematically by 



Fqq =TT + rS+ (55) 1 + {SS)2, 



(3) 



where T denotes thermal distribution and S shower dis- 
tribution. {SS)i signifies two shower partons in the 
same hard-parton jet, while {88)2 stands for two shower 
partons from two nearby jets, in a notation more self- 
explanatory than the one used in . For simplicity, we 
shall at times abbreviate (55) 1 by 55. 
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For pA collisions it may not be appropriate to refer 
to any partons as thermal in the sense of a hot plasma 
as in heavy-ion collisions. However, in order to main- 
tain the same notation for the decomposition in Eq. (jS)) 
independent of the collision types, we persist to use the 
symbol T to denote the soft parton distribution at low 
fcy, although they will occasionally be referred to as ther- 
mal partons, when it is more convenient. At low pT the 
observed pion distribution is exponential; we identify it 
with the contribution of the TT term. By writing T as 
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Cpi exp(-_pi/r). 



we obtain from Eq. Q 
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(4) 



(5) 



where T is the inverse slope. We shall determine C and 
T by fitting the d-|-Au data at low pT- 

For shower partons we can safely ignore the term {88)2 
in Eq. Q arising from two hard partons. S and {SS)i 
involve the partons of only one shower. They are convo- 
lutions of the hard parton distribution fi{k) with trans- 
verse momentum k and the shower parton distributions 
(SPD) Sj (z) from hard parton i to soft parton j (and j') 

m 

^(Pi) = E / dkkh{k)Sl{pi/k) , 

issupi,p2) = J2l dkkMk)i^Si[^),s{ 



P2 



k- pi 



The integrals begin at a minimum fcp below which the 
pQCD derivation of fi{k) is invalid. We set fco = 3 
GeV/c. The curly brackets in Eq. Q signify the sym- 
metrization of the leading parton momentum fraction 
[T^ . We have assumed in Eqs. and that the 
hard partons suffer no energy losses as they traverse the 
cold nucleus. 

The hard parton distributions fi{k) depend on the 
parton distribution functions (PDF) in the proton and 
nucleus, and on the hard scattering cross sections. For 
d-|-Au collisions Fries has performed the convolution and 
put the inclusive distribution for the production of a hard 
parton i at y = and at 0-20% centrality in a generic 
form over a wide range of k 



1 dcr''+'^" 
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where (Tin = 40.3 mb has been used. The parameters 
, fci , and Ui are given in Table I |3| • Nuclear shad- 
owing effects have been taken into account through the 
use of EKS98 PDF. The K factor is due to higher or- 
der corrections in pQCD. We shall set it at 2.5, as in 

illllii. 



TABLE I: Parameters in Eg. lISl fll. 
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We calculate the three contributions TT, T5, and SS 
to Fqq and then to dN^+^'^/pdp in Eq. I©. In the cal- 
culation there are two parameters: C and T. They are 
adjusted to fit the low pT region of the data. The point 
of view we adopt is that the soft component specified by 
C and T is not the predictable part of our model. In 
treating them as free parameters we do not compromise 
the predictable part of our model, which is the magni- 
tude of the contribution from the TS component in the 
recombination compared to the other components: TT 
at low pt and SS at high px ■ 

For d+Au collisions at RHIC, PHOBOS has published 
data on the pT spectra of charged particles, {h~^ + h~)/2, 
over the range 0.25 < pr < 6.0 GeV/c for various central- 
ities; however, the pseudorapidity range is 0.2 < 77 < 1.4 
[T^. PHENIX has preliminary data on 7r+ production 
over a narrower px range ( < 3 GeV/c) but for 77 = 
[l7l |. We use the latter because they are for 7r+ at 77 = 0. 
Although does not go above 3 GeV/c, the range is suf- 
(^cient to determine C and T, where the TT contribution 
dominates, and where the deviation from the exponential 
'jbehavior is just enough to reveal the TS contribution. 
{KDui prediction is the spectra for pr > 1 GeV/c. 
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FIG. 1: Pion distribution in transverse momentum compared 
to the data on 7r+ from PHENIX on d-|-Au collisions at 
-^5=200 GeV and 0-20% centrality. The three components 
in the recombination model are TT (light solid line), TS 
(dashed-dot line), and SS (dashed line). Heavy solid line is 
the sum of all three components. 

We show in Fig. 1 the three components together with 
their sum for 0-20% centrality. The agreement with the 
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data from ^3 is good, considering the fact that only 
the exponential component is adjusted to fit. The note- 
worthy features of our result is that the thermal-shower 
{TS) and shower-shower (SS) components both become 
more important than the thermal-thermal {TT) compo- 
nent above pr — 2 Gev/c and that TS is greater than 
SS for < 3 GeV/c. The values of the parameters 
determined are (for 0-20% centrality) 

C = 12 (GeV/c)"\ T = 0.21GeV/c. (9) 

With the above value of T we obtain from Eq. (O for 
the soft component of the pions (px = p) 

(pt) ^ 2T = 0.42 GeV/c, 

(pI) ^ 6T'^ = 0.26 (GeV/c)2, (10) 

while Eq. Q implies for the intrinsic transverse momen- 
tum kx of the partons (kx = Pi) 

(fcr) ^ T = 0.21 GeV/c, 

(4) = 2T'^ ^ 0.09 (GeV/c)2. (11) 

While the numbers in Eq. (|10|l are conventional, the in- 
trinsic width of the partons in Eq. Hll|) is very small 
compared to what is needed in the fragmentation models, 
generally (/c^) > 1 (GeV/c)^, even before broadening. 

The shower-shower component in Fig. 1 is dominant for 
Pt > 5 GeV/c; it is the same as the usual contribution 
from parton fragmentation [Til . The thermal-shower 
recombination is a unique feature of our model. It makes 
a dominant contribution in the 3 < pt < 8 GeV/c range 
in Au-I-Au collisions because of the large thermal com- 
ponent in the hot, dense system Here in the cold 
system only slightly excited, the values of C and T are 
lower, compared to 23.2 (GeV/c)-^ and 0.317 GeV/c, 
respectively, in Au-fAu collisions. Thus the TS contri- 
bution is subdued, but still large enough not only to cause 
a substantial deviation of the spectrum from exponential 
behavior, but also to give rise to the Cronin effect with- 
out large intrinsic (fc^) broadening, as we shall show. 
The ratio TS/SS is independent of the normalization of 
fi{k) and hence unaffected by the value of K. 

For other centralities of d-|-Au collisions we fix T at 
the value determined for 0-20% centrality, i.e., T = 0.21 
GeV/c, and use the values of (A^coii) given in 17] to 
rescale fi{k). We adjust the value of C to fit the low 
Pt normalization of the data [l7j . The results we obtain 
for dN^^^^ /pdp for all four centralities are shown in Fig. 
2 and are in good agreement with the data. The main 
point to stress is that the soft parton (fc^^ width remains 
the same at 0.09 (GeV/c)-^, without being broadened by 
successive kicks before hard scattering. In our approach 
what are different at more peripheral collisions are the 
decreasing values of C, which are 11, 7.8, and 5.65 for 
20-40%, 40-60% and 60-90%, respectively. The decrease 
of the density of the soft partons is a reasonable property 
of less central collisions, and can be generated in a Monte 
Carlo code; we merely take it from data, since it involves 




FIG. 2: Same as in Fig. 1 but for different centralities. The 
inverse slope T is the same in all cases; C is varied to fit the 
low-pT region. 

no new physics. The physics issue we want to emphasize 
is that the thermal-shower parton recombination is sen- 
sitive to the density of soft partons, and that component 
of the hadronization product affects the spectra in the 
moderately higher pT region, 1 < pT < 4 GeV/c. 

There are inaccuracies in our calculation due to the 
use of the lowest order pQCD results for the parameters 
in Table I, and the SPDs that inherit the uncertainties 
of the FFs Il3|. However, their effects tend to cancel, 
if we take the ratio of the calculated spectra at differ- 
ent centralities. PHENIX has preliminary data on the 
central-to-peripheral nuclear modification factor |l7j 

p , , (A^con)6o-90% dN^+^ypTdpTjO ~ 20%)^ ^^^ 

ncpiPT) = TTT ^12) 

(A^cou)o-2o% dNi^+^^/pTdpTi60 - 90%)' 

for Pt < Q GeV/c. We can determine Rcp{pt) directly 
from the results of our calculation; it is shown in Fig. 
3. Evidently, the theoretical curve agrees very well with 
the data [l3| . The agreement clearly lends support to our 
view that the Cronin effect is due to the recombination of 
soft and shower partons in the final state without initial- 
state broadening. 

The extension of our consideration to forward rapidi- 
ties naturally explains the decrease of the Cronin effect 
observed in d-|-Au collisions since the soft parton 
density decreases, as one moves from the central to the 
fragmentation region of the deuteron. We await the pre- 
cise data in the low-pT region at higher 77 in order to 
determine the 77 dependence of C in the soft component 
before calculating the pion spectra. 

On the basis that the thermal-shower recombination 
is the component mainly responsible for the enhance- 
ment of Kqp in the intermediate range of pt in Fig. 
3, we can infer qualitatively here that for proton pro- 
duction the corresponding BFqp should be even higher 
in the same region. The reason is that for three quarks 
to recombine in forming a proton, there are many more 
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FIG. 3: Central-to-peripheral nuclear modification factor. 
The data is from PHENIX for d-|-Au collisions at ^^=200 
GeV |l7|. The solid curve is the result of our calculation 
without adjustable parameter. 



terms: TTT + TTS + TSS + SSS + ■■■ [13, among 
which the TTS contribution has a quadratic depen- 
dence on C, and is therefore more sensitive to central- 
ity. Consequently, as the collisions become more central, 
dN^+^'^/pTdpT should receive a larger boost from TTS 
than does dN^'^^'^ /pxdpT from TS, resultin g in R^p 



being higher than i?Jp. The data of PHENIX [17| show 
that such a behavior has already been observed. That 
behavior is hard to interpret in a fragmentation model, 
since the broadening of the parton fc-p width in the ini- 
tial state should be independent of what a hard parton 
fragments into. To determine R^p quantitatively in our 



approach, we need to take the proton mass effects into 
account, for which our scale invariant formalism is inad- 
equate. That is a task to be done elsewhere. 

In conclusion we have shown, both here and in |l2j . 
that the separation of final states into independent and 
noninteracting soft and hard components is invalid, ex- 
cept when pt is very large. At moderate pt where the 
Cronin effect is found, the interaction between the soft 
and shower partons is important. Since the density of 
soft partons depends on the number of participants even 
in d-|-Au collisions, the hadron spectra at moderate px 
depend on centrality, when those soft partons recom- 
bine with the shower partons. Thus the enhancement 
of hadron production in more central collisions is a final- 
state effect, in contrast to the usual explanation in terms 
of initial-state fluctuations. The Cronin effect may now 
be regarded as another phenomenon in support of the re- 
combination model besides the large p/tt ratio and scal- 
ing elliptical flow 20]. Our result also gives credence to 
our assertion that shower partons form an essential com- 
ponent in the final state of a quark-gluon system pro- 
duced in a heavy-ion collision before hadronization takes 
place, but whose existence has hitherto been overlooked. 
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